Typical Scales in the Spatial Distribution of QSOs 



Zugan Deng^, Xiaoyang Xia^, Li-Zhi Fang^ 

Physics Department and Steward Observatory, University of Arizona, Tucson, Ari- 
zona 85721 

Graduate School, Chinese Academy of Sciences, Beijing, P.R.China 
Physics Department, Tianjing Normal College, Tianjing, P.R.China 



To appear in the Astrophysical Journal August 10, 1994 



Abstract 



We present results of searching for the possible typical scales in the spatial dis- 
tribution of QSOs. Our method is based on the second derivative of the two-point 
correlation function. This statistic is sensitive to the scale of the maximum in the spec- 
trum P{k) of the density perturbation in the universe. This maximum or bend scale 
can be detected as the wavelengths of the periodic component in the second derivative 
of the integral correlation function. For various QSO samples compiled from surveys 
of pencil-beam and bright QSOs, a typical scale of about 93 ± 10 h~^Mpc for go = 0.5 
has been detected. This typical scale is in good agreement with that found in the 
spatial distributions of galaxies, clusters of galaxies, and CIV absorption systems of 
QSOs if go is taken to be ~ 0.2. Therefore, it is likely a common or universal scale 
in the large scale structure traced by these objects. This result is consistent with the 
assumption that the typical scale comes from a characteristic scale in the spectrum 
of the density perturbation in the universe. 
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1. Introduction 



According to the standard scenario of structure formation in the universe, the 
initial perturbation produced by quantum fluctuation of scalar fields during the in- 
flationary era is scale- invariant. The power spectrum of the initial perturbation is 
assumed to be P{k) oc k"^, where k is the wavenumber of the perturbation, and the 
spectral index n ~ 1. Therefore, no typical scales exist in the very early universe. 
Subsequent evolution of the universe leads to a deviation of the density perturbation 
spectrum from a scale-invariant one. Typical scales emerge from the distribution of 
cosmic matter. For instance, in a linear regime, the density spectrum can be approx- 
imated as the following form (Peacock 1991; Mo et al 1993) 

where A is a typical scales, on which P{k) is the maximum, i.e. P{k) bends from oc k 
to at /c ~ 27r/A. The clustering of galaxies and clusters of galaxies showed that 
the bend scale A should be larger than about 100 h~^Mpc. On the other hand, the 
anisotropy of cosmic background radiation indicates that A < 1000 h^^Mpc. There- 
fore, in the standard model, the linear evolution of density perturbation brings out at 
least one typical scale in the range between 100 and 1000 h^^Mpc in the spectrum. 

The possible existence of typical scales in non-standard scenario of the structure 
formation has also been proposed. For a modified infiation model (Starobinsky, 1992), 
it has been found that a typical scale of about 100 h~^Mpc in the cluster- cluster 
correlation is crucial in determining the peculiarities of the inflation and the nature 
of the dark matter (Kotok et al 1993). 

Observations have indeed discovered structures in the distribution of galaxies or 
clusters with scales as large as about 100 h^^Mpc, including the great void (Kirshner 
et al. 1981), filaments and sheets (Haynes and Giovannelli, 1986), the Great Wall 
(de Lapparent et al. 1988), the Great Attractor (Dressier et al. 1987) and the 128 
Mpc 'periodicity' of pencil beam sample (Broadhurst et al. 1990). However, these 
observed scales cannot be identified as the bend scale in the density spectrum. 

In the last two years, systematic approaches to the typical scales in the spatial 
distribution of galaxies and clusters have been done by several groups. Buryak, et al. 
(1991, 1992) developed a method to probe typical scale from one-dimensional samples. 
Einasto and Gramann (1993) investigated the possible observational phenomena in the 
distribution of clusters and galaxies related to the bend scale in P{k). An extensive 
search for the typical scales has been made by Mo, et al. (1992a, b). Using the 
method of the second derivative of the integral two point correlation function, they 
have detected typical scales in the distribution of galaxies and clusters, especially, 
a scale of 130 h~^Mpc commonly exists in samples including the deep pencil-beam 
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survey (Broadhurst, et al. 1990), deep redshift surveys of Abell clusters (Huchra, et 
al. 1990) and QDOT survey of IRAS galaxies (Rowan- Robinson, et al. 1990). 

In this paper we extend this investigation to QSOs, i.e. searching for the possible 
typical scales in the spatial distribution of QSOs. Our motivations are twofold. First, 
the distributions of low-redshift objects like galaxies showed that the bend in the 
density spectrum may occur at wavelengths A ~ 150 h~^Mpc (Peacock 1991; Vogeley, 
et al. 1992; Peacock & West 1992; Mo, et al. 1992a,b; Jing & Valdarnini 1993). 
However, the wavelengths involved are already comparable to the sizes of the samples 
used, and the fair-sample assumption may then be questionable. This problem should 
be less severe for QSOs because we can have QSO samples with size much greater 
than 150 h~^Mpc. Secondly, if the typical scales detected in the structures of galaxy 
and cluster do come from the characteristic scales like the bend in the spectrum 
of perturbation, it should be measurable in the QSO distribution as well, unless 
QSOs trace substantially different large-scale structures than galaxies and clusters 
do. Therefore, it is important to see if the distribution of QSOs is consistent with the 
assumption that the typical scales found in galaxies and clusters are 'universal'. 

It is the interesting to note that the scale 100 h~^Mpc has already been mentioned 
in early studies of QSO clustering. About a decade ago, using the nearest neighbor 
analysis, Chu and Zhu (1983) showed that the distribution of QSOs hsted in the 
sample Bolton and Savege (1979) deviates from the Monte Carlo samples on the scale 
of about 100 h~^Mpc. Some authors also suggested the existence of isolated groups 
with comoving scales of about 100 h^^Mpc (e.g., Crampton, et al. 1989; Clowes & 
Campusano 1991). But these results do not provide a convincing argument for the 
scale considered. We will use more rigorous statistic to detect the typical scales in 
the samples of Boyle, et al. (1990, 1991) and Foltz, et al. (1987, 1989). 

Our plan is to give a brief description of the method in section 2, the results of 
typical scale analysis of QSO samples in section 3, a comparison of the typical scales 
of QSOs with that of galaxies in section 4, and a conclusion in section 5. 

2. Method of Detecting Typical Scales 

Statistics based on the amplitude of the two point correlation function ^(r) is the 
most popular method in the study of large scale structure. This method is, however, 
not adequate for detecting typical scales. The amplitude and the correlation length 
To in the two-point correlation function ^(r) = (r/ro)~^'^ do not relate to the bend 
scale in a simple way. For a given density spectrum P{k), the two point correlation 
function can be calculated by 



where B is a constant. Cenerally, the bend in the spectrum P{k) only leads to a 
slight drop in the amplitude of the correlation function on the scale of bending. This 




(2) 
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means, only the amplitudes of ^(r) at r ~ A are useful to probe A. However, on such 
large scales, the absolute value of the amplitudes of the two-point correlation function 
has a large statistical error due to the uncertainty in the mean density of objects 
considered. Structures on larger scales with density contrast less than the uncertainty 
of the mean density will be masked by the noise of the two-point correlation function. 
This problem is especially severe for QSO samples because the mean density of QSOs 
is redshift-dependent. Even with a homogeneous sample, it is still difficult to calculate 
the evolution of the mean volume density, because the deceleration parameter qo is 
poorly determined. 

The method developed by Mo, et al. (1992a,b) and Einasto and Gramann (1993) 
is based on the second derivative of the two-point correlation function. Wc will intro- 
duce this method by a slightly different way in order to demonstrate its advantage in 
searching for the typical scale, especially the bend scale in P{k). Let us consider the 
behavior of ^(r) when r is large. Eq.(2) shows that, for a spectrum with a maximum 
like that in eq.(l), the dominate term of ^(r) when r > A should be a periodic func- 
tion of r with wavelength equal to about A. For instance, if one takes an approximate 
form of the spectrum (1) as follows: P{k) = k for k < 27r/A and P{k) = k~^''^ for 
k > 27r/A, the dominant term of ^(r) at large r will be r"™ cos(27r/Ar), where m ~ 2. 
The second derivative of log,^(r) is then proportional to cos(27rr/A). Therefore, the 
bend scale A can be detected by the wavelengths of the periodic components in the 
second derivative (P log C,{r)/dr'^. Of course, such periodic components will also be 
masked by the noise given by the uncertainty of the mean density. However, it is well 
known from statistics that, for a given noise masked data set, identifying periodic 
components is easier than determining the absolute value of the amplitudes of the 
correlation function. Considering the sizes of QSO samples usually are much greater 
than the wavelengths of the periodic components involved, the statistic of detecting 
periodic component in a QSO sample would be more effective than that of determining 
the amplitude. 

In actual work, the usual two-point correlation function ^(r) is replaced by the 
function S(r) = 1 + ^{r), where ^(r) is the integral two-point correlation function 
defined by 

m = -s r^i^Wdx (3) 

When r > A, S(r) has about the same behavior as ^(r). Therefore, we can also use 
the statistic of (f logSi^r) / dr^ to detect the typical scales. In measuring clustering 
of high redshift objects on large scales, the statistic based on the integral two point 
correlation function ^(r) is sometimes more advantageous than ^(r). The reasons are 
as follows. 

First, for determining the two-point correlation function ^(r), one needs a choice 
of bin size of the separation of QSO pair. The binning may lead to false periodic 
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components in the ^ (r) with wavelengths equal to the harmonics of the bin scale, and 
then a misidentification of the typical scales. Moreover, the total number and number 
density in available samples of QSOs are very low, and the binning will cause a large 
fluctuation in ^(r) if the bin scale is chosen too small. These puzzles can be avoided 
by using the statistic ^(r) because, according to the definition eq.(3), it does not bin. 
Second, from eqs. (2) and (3), ^(r) can be related to the density spectrum by 



where the window function [sin kr — kr cos kr]/k dies off faster at large k than that for 
^(r) [eq.(2)]. Therefore, the statistical result will be less severely affected by clustering 
on small scales (Mo, ct al. 1993). 

When the boundary effect is negligible, 5(r) is given by 



where A^dd(r) is the number of QSO pairs with separation less than r, N^rir) is the 
mean number of object pairs between observed and random samples. A'" and Nr are 
the total numbers of objects in real and random samples, respectively. 

Eq.(5) shows that S(r) is given by an un- normalized integrated pair counts of 
QSOs, the result does not sensitively depend on the mean number density of QSOs. 
Therefore, the uncertainty in the mean density is avoided from the beginning. As a 
consequence, this method is not sensitive to the selection function used for generating 
the random samples as well. In this paper, we will fit one- dimensional samples by a 
cubic polynomial, and three-dimensional samples by a linear function in the redshift 
range of each sample. 

In calculating the derivative of S(r), we will meet differences like Ndd{T + Ar) — 
Ndd{i^)- Obviously, this difference will be dominated by noise when Ar is less than the 
mean distance D of nearest neighbor of QSOs in the sample. This will be the main 
source of the error in the derivative of the correlation function when Ar < D. In 
order to suppress the influence of this noise in small wavelengths, we smooth S(r) by 
convolution integral S(r) = JE{r')S{r — r')dr', where the smoothing function S{r) is 
equal to 1 when |r— r'| < L, and otherwise, and taking the smooth scale L to be equal 
to or less than D. Fluctuations with wavelengths less than the scale L will totally 
be suppressed in the function S(r) by the smoothing, while all inhomogeneities with 
scales comparable to or larger than the scale L will not be affected by the smoothing. 
Our algorithm is to use this smoothed function 5(r) to calculate the second derivative 
A6'(r) = (i^ log S(r)/(ir 2. 

The statistical significance of the peaks in the second derivative A^(r) can be 
measured by the standard deviation a which is estimated by Monte Carlo samples 




(4) 
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generated under the same selection conditions as the real samples. Usually we take 100 
random samples to calculate the standard deviation. Comparing the curve A6{r) of 
the real samples with that of random samples, we can infer the statistical significance 
of peaks appearing in the A6{r) of real samples. Estimating the significance in this 
way has the advantage that the edge effects are automatically avoided. 

The periodic components in A6{r) can be detected by power spectrum analysis 
(PSA). The wavelengths of these periodic components are the typical scales. The 
statistical significance of the existence of periodic components in the second derivative, 
A9{r), can be estimated by the usual way of power spectrum analysis. 

This method has been used to analyze 1- and 3-dimension samples of optical and 
IRAS galaxies or clusters of galaxies, and a common scale of 130 ±10 h~^Mpc was 
detected (Mo, et al. 1992a,b). The samples used for analysis include the deep pencil- 
beam surveys (Broadhurst et al. 1990), deep redshift surveys of Abell clusters (Huchra 
et al. 1990) and QDOT survey of IRAS galaxies (Rowan-Robinson et al. 1990). 
Therefore, the scale of 130 h~^Mpc might be a candidate for the bend scale A in the 
initial density spectrum P{k). Since the structures with scales as large as about 100 
h~^ Mpc in the present universe should still remain in the linear evolutionary stage, 
the typical scale found in the distribution of galaxies and clusters should probably 
also be measurable in the distributions of high redshift objects. Therefore, one should 
expect the existence of a 100 h~^Mpc typical scales in QSO distribution if QSOs trace 
the high peaks in the density field as galaxies and clusters of galaxies do. 

3. Statistical Results 

3.1 Pencil-beam samples 

The one-dimensional samples of QSOs used in our analysis are formed from Boyle 
et al 1990 (BFSP) and 1991 (BJS). BFSP contains about 420 QSOs identified in a 
complete (to B < 21), ultraviolet excess (UVX) survey, which covers 34 pencil beam 
fields, each has about 0.35 square degrees. These pencil-beam fields are scattered over 
eight 5° X 5° UK Schmidt fields. BJS includes 61 QSOs identified in a complete to 
B < 22 survey done by multicolored technique in three pencil-beam fields at high 
galactic latitudes. 

Redshift distribution of QSOs listed in BSFP and BJS are plotted in Figure la 
and b, respectively. As well known, the UVX and multicolored technique are likely 
to provide QSO candidates with high completeness when redshift z < 2.2 (Veron 
1983). On the other hand, the imposed stellar morphological criterion may cause 
the incompleteness at redshift less than about 0.6. Figure 1 shows that most QSOs 
in the BFSP and BJS are in the redshift range of 0.6 to 2.2 and the number of 
QSOs dramatically decreases outside this interval. Therefore, the samples consisting 
of QSOs with redshifts from 0.6 to 2.2 in BSFP and BJS should be largely complete 
and unbiased. We adopt, respectively, QSOs with 0.6 < z < 2.2 in 1) BSFP and 2) 
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BJS as two parent samples in our statistic. 

Because each field covers only about 0.35 deg^, the size of their cross section is 
about 20 h^^Mpc at z 1. One can consider these samples as one-dimensional if we 
focus on the structures with scales much larger than 20 h~^Mpc. Each one-dimensional 
sample can be seen as a representation of the three-dimensional distribution in a given 
direction. Obviously, some features shown in these samples arc direction-dependent. 
In order to reduce the influence of the local features on the statistics, we use combined 
sample, which consists of a number of the pencil-beam samples in different directions. 
For such combined samples, the individual features of the pencil-beam fields should 
be less important. 

We made four subsamplcs called QN, QS, QSGP and QBJS, which consist of 
QSOs in the fields of northern sky, southern sky (excepted those in the Schmidt field 
along southern galactic pole direction), southern galactic pole and that given by BJS 
survey, respectively. Table 1 shows the numbers of pencil-beam fields Np, numbers of 
QSOs Nq for each subsample. It should be pointed out that a Schmidt field is about 
5° X 5°, two pencil-beam fields in the same Schmidt plate have a mean separation 
of about 2°. 5 which corresponds to an across scale of about 90 h^^Mpc at 2; ~ 1. 
Thus, in studying the structures with scales larger than 100 h~^Mpc, the pencil-beam 
samples in the same Schmidt plate should not be considered as totally independent 
samples, i.e. different pencil-beams may imprinted by a same structure on scale larger 
than the separation of the pencil-beams. 

Figures 2a-d present the results of the second derivative, A^^(r), for each subsam- 
ple. The bold lines show the A9{r) for the real samples, and the light lines are the 
standard deviations ilcr given by random samples. Since the mean distance of near- 
est neighbor QSOs in these samples is equal to or larger than about 50 h~^ Mpc, the 
smoothing scale L is taken to be 40 h~^Mpc. The sharp dips appearing in Figure 2 on 
scales less than about 40 h~^Mpc are caused by that, for all random samples, N^rir) 
is zero on small scales. It can be clearly seen from Figure 2 that, for all subsamples of 
QSOs, the curves of A^(r) show periodically distributed peaks (and valleys). Figure 
3a and b show the same calculation as Figure 2 but the smoothing scale L is taken 
to be 20 h^^Mpc. As expected, more fiuctuations with short wavelengths appeared 
in Figure 3, while the fiuctuations on scales larger than 40 h~^Mpc are the same as 
the case of L = 40 h~^Mpc. 

The statistical significances of each peak (and valley) in A^(r) are marginally 
higher than ±0", mostly in the range of 1 — 2.5 a. This result is the same as that 
given by previous studies. Up to now, almost all QSO structures detected with scales 
greater than 10 h~^Mpc are in the significance level of 2 - 3 cr (Ivovino & Shaver, 
1988; Bahcall & Chokshi 1991; Boyle & Mo 1993). However, our method is not only 
based on the individual peak in A^(r), but in the regular or periodic distribution 
of these peaks (and valleys). Figure 4 and 5 plot, respectively, the power spectrum 
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of the Ae{r) for sample QS and QBJS with L = 40 and 20 h'^Mpc. One can find 
from these figures that 1) all power spectrum show a peak around A ~ 90 h~^Mpc 
with confidence no less than 99% (for L= 40 h~^Mpc); 2) if we use the width of the 
peak as a measure of the uncertainty of the wavelength, the mean wavelength for 
one-dimensional samples is 90 ± 8 h~^Mpc; 3) for each subsample the spectrum of 
L = 20 and 40 h~^Mpc have the same shape, therefore these statistical results are 
independent of the smoothing length. 

3.2 Three dimensional samples 

The three-dimensional sample used in this paper is compiled from the LBQS 
survey (Foltz et al., 1987, 1989; Hcwctt ct al., 1991; Chaffee et al., 1991; Morris ct 
al., 1991). The LBQS survey presented more than 1000 QSOs with rrij <18.5. The 
total area of these fields is about 800 square degrees. It is one of the largest and 
uniformly selected QSO sample up to date. The redshift interval of this sample is 
between 0.2 and 3.3. An artificial cut-off has been made at redshift 0.2, because too 
many stars mixed into the candidates below this redshift. As the rcdshifts of QSOs 
becomes higher than 3.3, the Lya line will move out of the j band. The sample can 
be considered to be complete in the magnitude interval 16.0 < rrij < 18.7. Figure 6a 
plots the redshift distributions of QSOs in samples LBQS. 

The first sample compiled from this survey is called LBQS, which contains of all 
LBQS QSOs with redshift in the range from 1.0 to 2.2. The limitation of 1.0 < z < 2.2 
comes from the following consideration, a) Each plate used in LBQS survey coves an 
area of about 6° x 6°, which spans ~ 200 h~^Mpc and higher when z > 1 and go = 0.5. 
Therefore, if we are interesting in probing structures with scales equal to or larger than 
100 h^^Mpc, only the sub-samples with z > 1 can be treated as a three dimension one. 
b) The number of the LBQS QSOs drops rapidly when z > 2.2. The mean distance of 
nearest neighbor QSOs at 2; > 2.2 is equal to or larger than ~ 100 h~^Mpc. Therefore, 
the data at z > 2.2 are no longer suitable for probing structures with scales of about 
100 h-^Mpc. 

To study the possible infiuence of the foreground objects on the typical scales, 
such as that given by gravitational lensing effect, we compiled a sample called LBQS- 
V, which consist of all LBQS QSOs with 1.0 < 2; < 2.2 excepting those in field of the 
nearest supercluster Virgo (Figure 6b). The redshift distributions of both LBQS and 
LBQS-V are quite smooth. Samples LBQS and LBQS-V are listed in Table 2. 

Figures 7a and b plotted the result of A6{r) for the LBQS and LBQS-V, respec- 
tively. The ranges of ilu given by random samples are shown by the light curves. 
The ratio of signal to noise in 3-dimensional samples appears to be higher than that 
of the pencil beam samples. The significance of the peaks (and valleys) now is about 
3(7. As in the 1-dimensional sample, the significant peaks distributed regularly or pe- 
riodically. Figures 8 and 9 show, respectively, the power spectrum of samples LBQS 
and LBQS-V with L = 40 and 20 h^^Mpc. The mean wavelength is ~ 95 ±9 h^^Mpc. 
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This is the same as that of pencil beam samples. 

The value of A ~ 90-100 h~^Mpc found here is in agreement with those obtained 
from the amplitudes of the correlation function (Mo & Fang 1993). Using a best 
fit of the integral correlation function to the power spectrum P{k) [eq.(l)], it found 
A ~ 100 — 200 h~^Mpc. Therefore, the scale of ~ 100 h~^Mpc seems to be universal 
for the various QSO samples considered. 

4. Difference of Typical Scales between QSOs and Galaxies 

The typical scale found in QSO distribution shows a difference from that of galax- 
ies and clusters of galaxies (Mo, et al. 1992a), which was found to be 130 ±10 h~^Mpc 
by the same method. It is generally believed that the structures with scales larger 
than about 50 h~^Mpc should still remain in linear evolution regime. If the typical 
scale comes from the bend scale in the initial density spectrum, the comoving typical 
scale of QSOs should be the same as that of galaxies and clusters of galaxies. There- 
fore, it is necessary to study the possible origin of the observed difference between the 
typical scales of QSOs, and galaxies and clusters. 

If the typical scale is assumed to be local, i.e. being only a feature of nearby 
galaxies, one should not expect that the same typical scale shows up in the spatial 
distribution of QSOs. However, it has been found that the typical scale of 130 h~^Mpc 
exist in almost all samples of galaxies and clusters of galaxies (Mo, et al, 1992a, b). 
On the other hand, the fact that no difference has been found between the results of 
LBQS and LBQS-V (Figure 6a and b) indicates that the 95 h~^Mpc typical scale of 
QSOs may not be affected (at least in current error bar) by gravitation lensing of a 
local structures like the Virgo clusters. Therefore, one cannot explain the difference 
of the typical scale between galaxies and QSOs as a local effect. 

If galaxies and QSOs trace different aspects of the density fields in the universe, we 
should not expect that galaxies and QSOs have the same typical scales. Considering 
the bias mechanism for galaxies is probably no longer useful for QSO formation, 
one may reasonably assume that the QSO-traccd structures are different from that 
traced by galaxy. Indeed, the number density of QSOs is much less than galaxies. 
Therefore, in terms of bias model, the biasing threshold of QSO should be higher than 
galaxies. Bower, et al. (1993) recently proposed a bias model describing cooperative 
formation of galaxies, in which the threshold of galaxy formation is scale-dependent, it 
is lower in a domain with higher mean density, and higher in the area with lower mean 
density. This is equal to replacing P{k) by P{k)B{k), where the bias function B{k) 
is decreasing with k increasing. Obviously, the bend scale of "spectrum" P{k)B{k) 
will be greater than that of P{k). 

However, this explanation encounters difficult if we consider the following facts. 
First, QSO clustering satisfies the same power law correlation function as galaxies 
and clusters (Shanks et al. 1988; Fang ct al. 1985; Shaver 1988; Chu & Zhu 1989; 
Crampton, Cowley & Hartwick 1989; Boyle 1991). Second, QSO clustering is the same 
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as small groups of galaxies or poor clusters (Bahcall & Chokshi 1991). It has been 
known for a decade that low redshift QSOs are preferentially located in poor clusters 
or groups of galaxies. This was found by the QSO-galaxy covariance function (Yee & 
Green 1987), CIV-associated absorption in high redshift radio-loud QSOs (Flotz et al 
1988), clustering analyses of the QSO distribution and galaxy environments around 
of QSOs (EUingson ct al 1991a). It has also been shown that the velocity dispersion 
of galaxies around QSOs is ~ 400 km s^^ (EUingson et al 1991b). This means, QSOs 
trace the same density field as poor clusters do. Therefore, the formation and radiation 
of QSOs may not provide an effective mechanism that leads to the difference of typical 
scales between QSOs and clusters of galaxies. 

Now we turn to the explanation based on go-dependence of the typical scale. In 
the previous sections, all scales are calculated under the assumption that the universe 
is of Einstein-de Sitter, and thus the deceleration parameter qo is taken to be 0.5. As 
it has been pointed out by Shank et al (1987) and Mo, et al. (1992a, b), the typical 
scale is crucially dependent on qo for high redshift objects like QSOs. Figures 10 is the 
power spectrum of the A9{r) for QBJS when go ~ 0.2, 0.4 and 0.7, respectively. From 
the peaks in Figure 10 one can see a systematic increase of the wavelength with the 
decrease of go- This relationship is plotted in Figure 11. It shows that when go = 0.2, 
the typical scale of QSOs is 125 ± 11 h~^Mpc. For other samples, we found about 
the same result. Figure 12 and 13 showed that all power spectrum of samples QN, 
QSGP, LBQS and LBQS-V have a peak at 130 ± 15 h^^Mpc. 

Since galaxies and clusters have low redshifts, their typical scales do not depend 
on go. Therefore, the typical scale of QSOs is in good agreement with that of galaxies 
and clusters (Mo, et al. 1992a, b) when go ~ 0.2. In other words, if one assumes that 
the comoving value of the typical scale can be used as a "standard cosmological rod" 
in the linear regime of an expanding universe, the universe should be of go ~ 0.5. 
Considering various uncertainties in the typical scales of galaxies, clusters, and QSOs, 
it would be better to say that go should be less than 0.5 at 95% confidence. 

It is interested to point out that the conclusion of a low total mass density universe 
(go ~ 0.1) has also been proposed by several independent researchers (Park, et al. 
1992; Bahcall and Cen, 1992; Vogeley, et al., 1992). In the paper by Mo, et al 
(1992b), the go-dependence of the typical scales of the CIV absorption systems and 
the Lya forests of high redshift QSOs has been studied. Using their result, one can 
find that the value of go, at which the typical scale of CIV system is the same as that 
of galaxies and cluster, is go ~ 0.2 — 0.3. Therefore, when go = 0.2, the 130 h~^Mpc 
typical scale is likely universal among the samples of galaxies, clusters, QSOs and 
the CIV absorption systems of QSOs. However, for this go, the typical scale of Lya 
forest is different from 130 h~^Mpc. This is understandable because the CIV systems 
probably originated from absorption of clouds associated with galaxies (Waymann, et 
al. 1979; Young et al, 1982), and Lya comes from clouds which are unable to form 
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QSOs and galaxies. In a word, it would be reasonable to say that the 130 h ^Mpc 
typical scale is universal in the distributions of galaxies, clusters and quasars. 

5. Conclusion and Discussion 

The possible typical scales in the distribution of QSOs have been detected by 
means of the second derivative of integrated correlation function. A typical scale of 
about 93 ± 10 h~^Mpc when go = 0.5 have been detected with considerable confidence 
in available 1- and 3-dimension samples of QSOs. This typical scale is probably 
"universal" for various subsamples of QSOs. If qo is taken to be 0.2, the QSO typical 
scale becomes the same as that of galaxies and clusters of galaxies. One can then 
have the following conclusions. 

1. The existence of a ~ 100 h~^Mpc typical scale in the spatial distribution of 
QSOs further strengthens the picture that QSOs probably trace the same larger scale 
density field of the universe as galaxies and clusters of galaxies do. 

2. The redshifts range and sizes of QSO samples used here are totally different 
from that of galaxies and clusters of galaxies. The agreement of the QSO typical scale 
with the result of galaxies and clusters suggests that the detected typical scale is most 
likely universal in the large scale structure. This is consistent with the assumption that 
the typical scale is due to a characteristic scale in the initial perturbation spectrum, 
such as the bend scale in the density spectrum at the linear regime. 

3. If the typical scale comes from the initial spectrum of the density perturbation, 
one can use the comoving values of the typical scale as a cosmic standard length when 
this scale remains in linear evolution. Accordingly, if we require that the comoving 
value of typical scale of QSOs is the same as that of galaxies, clusters and CIV 
absorption systems, the typical scale measurement favors an open universe, i.e. go < 
0.5. 

4. The bend scale is model-dependent. For instance, CDM model should have a 
lower bending scale than that of a hybrid model. The formation time of structures 
with scale as large as the bend scale is also model-dependent. Therefore, the fact that 
the structures of bend scale exist at the high redshift seen for QSOs distribution may 
help discriminate among various cosmological models. 
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useful discussion. Deng is grateful to the University of Arizona for hospitality and 
support from the international exchange program. D and X thank also the support of 
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Table 1 

Data of 1-dimensional samples 



Sample Np Nq Notes 



QN 
QS 



13 
14 



141 northern sky (BSFP) 
154 southern sky excepting 



QSGP 
QBJS 



7 
3 



southern galactic pole (BSFP) 
94 southern galactic pole (BSFP) 
52 (BJS) 



Tabic 2 



Data of 3-dimensional samples 



Sample Np Ng Notes 



LBQS 18 510 Foltz et al, 1987, 1989; Hewett et al. 1991 and 

Chaffee et al., 1991, Morris et al. 1991 

LBQS-V 14 399 excepting Virgo fields 

(references are the same as LBQS) 
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Figure captions: 

Figure 1 Redshift histograms of QSOs listed in pencil beam surveys, a. BFSP given 
by Boyle et al. (1990), and b. BJS given by Boyle et al. (1991). 

Figure 2 Curves of A^(r) of samples a) QN; b) QBJS; c) QS; and d) QSGP. Light 
lines show the curves of ±0" given by a average of 100 random samples, qo is 
taken to be 0.5, and smoothing scale L — AO h~^Mpc. 

Figure 3 Curves of A9{r) of samples a) QN and b) QS. These curves are obtained 
in the same way as Figure 2, but taking the smoothing scale L — 20 h~^Mpc. 

Figure 4 Power spectrum of A6{r) for sample QS. The smoothing length L is taken 
to be 40 and 20 h~^Mpc, and qo — 0.5. For peaks with P > 7.6, the confidence 
of the existence of a periodic component is > 99%. 

Figure 5 Power spectrum of A9{r) for sample QBJS. 

Figure 6 Redshift distributions of QSOs hsted in the LBQS survey (Foltz et al. 1988, 
1989): a. LBQS, consisting of all LBQS QSOs; b. LBQS-V, consisting of all 
LBQS QSOs, but excepting those in the area of the Virgo cluster. 

Figure 7 Curves of A6{r) of three-dimensional samples, a. LBQS; b. LBQS-V, 
respectively, qo is taken to be 0.5. 

Figure 8 Power spectrum of A9{r) for sample LBQS. 

Figure 9 Power spectrum of A9{r) for sample LBQS-V. 

Figure 10 Power spectrum of A6{r) for sample QBJS. The deceleration parameter 
is taken to be qo — 0.2, 0.4, 0.7, respectively. 

Figure 11 QSO's Typical scale in sample QBJS as a function of go- 
Figure 12 Power spectrum of A9{r) for samples QN and QSGP when qo — 0.2 and 
L = 40 h-^Mpc. 

Figure 13 Power spectrum of A6{r) for samples LBQS and LBQS-V when qo — 0.2 
and L = 40 h'^Mpc. 
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